Upper Bound of 0.28eV on the Neutrino Masses from the Largest Photometric Redshift Survey 
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We present a new upper limit of ^ m v < 0.28 (95% CL) on the sum of the neutrino masses assuming a flat 
ACDM cosmology. This relaxes slightly to J2 m " < 0-34 and J2 m " < 0-47 when quasi non-linear scales 
are removed and w ^ — 1, respectively. These bounds are derived from a new photometric redshift catalogue 
of over 700,000 Luminous Red Galaxies (MegaZ DR7) with a volume of 3.3 (Gpc /i -1 ) 3 , extending over the 
redshift range 0.45 < z < 0.65 and up to angular scales of ^ max = 300. The data are combined with WMAP 
5-year CMB fluctuations, Baryon Acoustic Oscillations (BAO), type la Supernovae (SNe) and an HST prior on 
the Hubble parameter. This is the first combined constraint from a photometric redshift catalogue with other 
cosmological probes. When combined with WMAP this data set proves to be as constraining as the addition 
of all SNe and BAO data available to date. The upper limit is one of the tightest and 'cleanest' constraints 
on the neutrino mass from cosmology or particle physics. Furthermore, if the aforementioned bounds hold, 
they all predict that current- to-next generation neutrino experiments, such as KATRIN, are unlikely to obtain a 
detection. 
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Introduction - Studies of the neutrino have traditionally 
been the realm of particle physics experiments, with Super- 
Kamiokande Q first indicating the presence of mass. Neu- 
trinos were shown to oscillate between the known flavors (y e , 
v^, v T ) solving, in the process, the long standing solar neu- 
trino problem. This implies the neutrinos have at least two 
non-zero mass eigenstates (mi, 777,2, 777-3) because the fla- 
vor mixing depends on the differences between their masses 
squared. Subsequently, bounds have been placed on the split- 
ting between the neutrino mass eigenstates from a host of 
solar, accelerator and atmospheric experiments; | A 777,3 -J ~ 
2.4 x 10- 3 eV 2 and Ara^ « 7.7 x l(T 5 eV 2 (e.g. Q). How- 
ever, currently both the absolute scale and the hierarchy of 
the masses remain hidden. KATRIN, a kinematic beta decay 
experiment [3 ], aims to provide a constraint in the future. 

Cosmology not only probes the absolute mass scale of 
the neutrino but is a completely independent method to test 
against e.g. (H EL In any case, it is imperitive to include 
an accurate prescription for the neutrino in cosmology, as any 
failure to do so can bias the other cosmological parameters. 
A cosmological constraint on the sum of the neutrino masses 
is primarily a constraint on the relic Big-Bang neutrino den- 
sity Vt v . One can relate this density to the sum of the mass 
eigenstates m v (e.g. (6)) as given by, 

n - a) 

"~ 93.14h 2 eV Uj 

The direct effects of the neutrinos depend on whether they are 
relativistic, non-relativistic and the scale under consideration. 
Neutrinos have a large thermal velocity as a result of their 
low mass and subsequently erase their own perturbations on 
scales smaller than the free streaming length (5J0. This sub- 
sequently contributes to a suppression of the statistical clus- 
tering of galaxies over small scales and can be observed in a 
galaxy survey. The abundance of neutrinos in the Universe 
can also have a direct effect on the primary CMB anisotropics 
if non-relativistic before the time of decoupling (i.e. when 
sufficiently massive). However, one of the most clear effects 



at this epoch is a displacement in the time of matter-radiation 
equality. All these cosmological effects can be used to impose 
bounds on the neutrino mass. 

Previous studies have capitalised on these signatures and 
have started to place sub-eV constraints on the absolute mass 
scale l8lfT3l We utilise the new SDSS MegaZ LRG DR7 
galaxy clustering data lITH to provide the first photometric 
galaxy clustering constraint on the neutrino and, combining 
with the CMB, examine the complementarity of these early 
and late-time probes. With an almost comprehensive combi- 
nation of probes this renders one of the tightest constraints on 
the neutrinos in cosmology and therefore physics. 

Assumptions - We assume a flat Universe with Gaussian 
and adiabatic primordial fluctuations and a constant spectral 
index. The effective number of neutrinos are fixed to N e ^ = 
3.04, thereby assuming no sterile neutrinos or other relativistic 
degrees of freedom. The constant dark energy equation of 
state is at first set to w = — 1 and later relaxed. Finally, we 
consider the neutrinos to be completely mass degenerate given 
that current inferred bounds are much greater than the splitting 
hierarchies. The potential of future surveys to discriminate the 
mass hierarchy has been discussed in fT5U20ll . 

Analysis - Although parameter degeneracies and a mild 
insensitivity to relativistic (lighter) neutrinos limit the upper 
bound one can place on m v w ^ m m e CMB l2T1l it repre- 
sents a clean and relatively systematic-less cosmological tool 
whose high statistical discrimation of the remaining cosmo- 
logical model facilitates a competitive combination of probes. 
We therefore start by using the latest 5 -year WMAP data and 
likelihood as described in l22l to vary seven ACDM parame- 
ters: Q b h 2 , Q c h 2 , ^a, n s , r, ln(10 10 A s ) and Asz, in addition 
to m u . r, n s and A s represent the optical depth to reionisa- 
tion, the scalar spectral index and the amplitude of curvature 
pertubations defined at k = 0.002/Mpc, respectively. The 
contributions from the Sunyaev-Zeldovich fluctuations are in- 
cluded by adding a template spectrum Cf z with pre-f actor 
Asz following [23]. This is allowed to vary as < Asz < 2 
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FIG. 1 : The best fit angular power spectra Ci in the combined anal- 
ysis (solid lines) are plotted over the MegaZ DR7 data. The panels 
relate to four redshift bins with width Az — 0.05 from z — 0.45 
(main panel) to z = 0.65 (panel 4). These spectra are good fits 
to the galaxy statistics including scales not utilised in the analysis 
(£ > 300). The fit is also plotted for linear spectra (dashed lines) and 
highlights the scale at which the non-linear regime starts to become 
significant. The dotted line demonstrates the effect of introducing 
m v = 1 eV neutrinos with all parameters, except Q c , held fixed. 

1221 . We use the pre-March 2008 version of CAMB (21 to 
produce the CMB power spectra. The reionisation is there- 
fore treated as a semi-instantaneous process. The gravitational 
lensing effect on the CMB is also included, e.g. l25ll , and the 
CosmoMC package l26l is used for parameter exploration. 

Our CMB run yields ]T m v < 1.271 eV at the 95% con- 
fidence level consistent with Komatsu et al. fT2l . This bound 
implies the neutrinos were relativistic at decoupling and as 
such induces a degeneracy between the neutrino masses and 
ft m as well as the Hubble parameter h. This can be seen in 
Figure [2] as well as fi"fl[T2ll2ll . This degeneracy can be im- 
proved by adding supernovae (SNe) data from the first year 
Supernova Legacy Survey (SNLS (27)) and the BAO data 
from [28]. Our analysis for WMAP + SNe + BAO gives 
Y, m v < 0.695 eV (95% CL) similar to Komatsu et al. C3 
(E m v < 0.67 eV) and Ichiki et al. El (£ rn v < 0.76 eV). 
In order to go beyond such studies we include the MegaZ LRG 
(DR7) photometric redshift survey that will be presented in 
lfT4l . which we have checked to be compatible with earlier 
SDSS clustering (22 andphoto-z analyses (30l[3T]|. This adds 
galaxy clustering information that is sensitive to the growth 
of structure suppressed by the free streaming neutrinos. The 
SDSS colours provide reliable photometric redshift estimates 
and, due to their high luminosity, probe a large region of cos- 
mic volume. Encapsulating 7746 deg 2 , we utilise 723,556 



photometrically determined LRGs in four redshift bins of 
width Az = 0.05 between 0.45 < z < 0.65 in a spherical 
harmonic analysis of the galaxy distribution until a maximum 
multipole I max — 300. These galaxies are calibrated by the 
2SLAQ redshift survey E2 using ANNz (S3 as in ED and 
the previous DR4 photometric release (29l|33]|. Specifically 
we use the angular power spectrum defined as, 

d =< 5 2D 5* 2D >= 4tt J A 2 (k)W 2 (k)^. (2) 

where A 2 (k) is the dimensionless power spectrum calculated 
with CAMB. The matter distribution is projected onto a plane 
in the sky with weight W 2 (k) in this statistic described by 
both Wi(k) = J f{z)ji{kz)dz and f(z) = n(z)D(z)(^) 9 
with the spherical Bessel function ji(kz), the linear growth 
factor D(z) and the normalised redshift distribution n(z). The 
effects of redshift space distortions are included as described 
in (34l[35l. The likelihood combines the four measured red- 
shift bins and includes the full covariance as a result of pho- 
tometric errors scattering galaxies between bins and therefore 
correlating slices. There are four additional parameters in- 
cluded in the study as a result of the galaxy bias in each bin 
0>i, 62, 63 and 64), i.e. modestly accounting for the redshift 
dependence, which is seen to increase at high z (141 |29) . De- 
spite the non-linear contribution becoming significant only at 
scales t > 300 we use halofit [36] to model the non-linear 
power spectrum. It is interesting to note that the point corre- 
sponding to the largest angular scale band in the highest red- 
shift bin indicates an excess of power. Hints of this were seen 
in the earlier photometric releases by (221 E) with the excess 
labelled by the arrow in panel 4 in Figure [T] 

This survey is not only one of the most recent and largest 
to date but is one of the most competitive available. How- 
ever, these power spectra provide an additional incentive for 
this combined measurement. This because the BAOs, which 
were shown to be so advantageous before, can be used in con- 
junction to MegaZ with no cross-covariance. The BAO data is 
extracted at z = 0.2 and z = 0.35, whereas MegaZ is defined 
at a higher redshift. They therefore constitute two independent 
data sets and can be used both simply and simultaneously. 

By combining the MegaZ LRGs as described above with 
the previous CMB, SNe and BAO data in a complete joint 
analysis we find a significantly lower bound of ^2m v < 
0.325 eV at the 95% confidence level. Again, this is roughly a 
factor 2 improvement in the neutrino masses with the addition 
of the LRGs and is shown clearly against Q m , h and the ID 
marginalised distribution in Figure [2] 

The information on the growth of stucture is paramount to 
the improvement seen in this study. However, part of this 
information originates from the quasi-non-linear regime and 
could systematically bias the inferred constraint. While ap- 
proaches are developed and work continues into the effects of 
the neutrino on these scales (E.g. l37l ) we repeat the com- 
bined analysis with the smaller scales removed. By trun- 
cating the multipoles at ^ max = 200 this more conservative 
approach is seen to give a similar but slightly relaxed limit 
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FIG. 2: Left Panel: The marginalised ID distribution for the neutrino from three incrementally combined analyses. The vertical dashed 
lines correspond to 95% confidence levels. Other Panels: 68% and 95% marginalised distributions for the matter density Q m , Hubble 
parameter h and spectral index n s against the total neutrino masses. The contours correspond to (from bottom layer) WMAP-only (red/dark), 
WMAP+MegaZ (yellow/light), WMAP+SNe+BAO (blue/dark) and WMAP+SNe+BAO+MegaZ+HST (green/light). Table |l| highlights that 
the vast majority of the gain in the last two contours originates from the new and complementary galaxy clustering data. 



of m v < 0.393 eV. While this highlights the importance 
of understanding non-linearities for obtaining the most strin- 
gent constraints, it is reassuring that there is still a marked 
improvement on the previous study (CMB+SNe+BAO) with 
linear LRGs. 

It is also intriguing to compare the input of the LRGs to 
those of the two distance measures (SNe+BAO). These have 
previously been highly beneficial to the uncertainty. We there- 
fore perform a joint analysis using just the WMAP5 and LRG 
data, subsequently obtaining the limit ^m v < 0.651 eV at 
the 95% confidence level. This is comparable to the spec- 
troscopic DR7 galaxy clustering addition to the CMB in 031 
with ^ m v < 0.62 eV and illustrates the development of pho- 
tometric surveys as a competitive tool for the future. 

We conclude this work by further restricting the cosmolog- 
ical parameter space with the addition of the new HST prior 
on the Hubble parameter to the WMAP5 + SNe + BAO + 
MegaZ DR7 run. The improved prior was recently found to 
be: H = 74.2 ± 3.6 km s^Mpc" 1 [ 38 1. With this the final 
limit in this study is reduced to ^ m v < 0.28 eV at the 95% 
confidence level. This is one of the tightest constraints in the 
literature. The angular power spectra Ci corresponding to the 
best fit values are plotted in Figure [T] with the galaxy cluster- 
ing data. An overview of all the neutrino bounds are displayed 
in Tableland a plot of all parameter combinations compared 
to the CMB -only study is displayed in Figure|3] Our estimates 
on the bias are h = 1.74 ±0.07, b 2 = 2.02 ±0.08, b 3 = 2.12 
±0.09 and b A = 2.39 ± 0.10. 

For w ^ — 1 the tighter bound relaxes slightly to Y m v < 
0.47 eV, which should be compared to m v < 0.62 eV from 
l39l . We note that biasing could act to mimic the neutrino sig- 
nature over smaller scale analyses. As a gauge of this effect 
we implement, as an example, the 'Q-model' of l40l ; result- 
ing in a combined constraint (all data) of Y m v < 0.44 eV. 
However, undertaking the challenge of modelling the possi- 
bility of scale dependent bias accurately and self consistently 
is left to future work. 



E m u (95% CL) Analysis 



< 1.271 eV WMAP5 

< 0.695 eV WMAP5 + SNe + BAO 

< 0.651 eV WMAP5 + MegaZ 

< 0.393 eV WMAP5 + SNe + BAO + MegaZ (€20o) 

< 0.344 eV WMAP5 + SNe + BAO + MegaZ ( ^ 20o) + HST 

< 0.325 eV WMAP5 + SNe + BAO + MegaZ 

< 0.281 eV WMAP5 + SNe + BAO + MegaZ + HST 

< 0.491 eV WMAP5 + SNe + BAO + MegaZ (€20o) + HST 

< 0.471 eV WMAP5 + SNe + BAO + MegaZ + HST 

TABLE I: A summary of the bounds placed on ^ tti u in this letter. 
£200 corresponds to the truncation in the maximum multipole scale 
to remove the quasi non-linear regime. The top constraints are for 
w — —1; the bottom for iu / —1, marginalised over. 



Conclusions - Using the biggest ever large scale structure 
survey we have set bounds on the neutrino masses at m v < 
0.28 eV (£ max = 300) and^m, < 0.34 eV (£ max = 200) 
at 95% CL, when combined with WMAP5+SNe+BAO+HST 
data. This is the first ever determination of neutrino masses 
from a photometric galaxy redshift survey. Not only have we 
shown that photometric redshifts can be used for this prob- 
lem, but also that such a galaxy survey is competitive with all 
currently available geometric probes (SNe+BAO) or spectro- 
scopic clustering when added to the CMB. Our constraint is 
one of the tightest current bounds available without the use 
of data from Lyman-a (e.g. |9 |), which is prone to systemat- 
ics, or a complicated modelling of the bias BH . Further, all 
our results show that KATRIN's 1 3 ] projected 90% sensitivity 
m v ^0.6 e V) leaves an unlikely neutrino mass detection. 
Finally, our overall method shows great promise for the next 
generation of surveys, which will yield upper limits of e.g. 
0.12 eV El and 0.025 eV H2 at 95% CL. 
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FIG. 3: The 2D 68% and 95% contours and marginalised ID dis- 
tributions for 7 cosmological parameters (Qbh 2 , Q c h 2 , Q\, n s , r, 
ln(10 10 A s ) and £ m u ) in a WMAP5 + SNe + BAO + MegaZ DR7 
+ HST combined constraint (red contours). The amplitude of the 
Sunyaev-Zeldovich fluctuations (Asz) and four bias parameters (&i, 
&2, &3, &4) are marginalised over but not plotted. The black contours 
are given for the WMAP-only analysis. 
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